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ABSTRACT 

We report on the long term X-ray monitoring with Swift, RXTE, Suzaku, Chandra, and XMM- 
Newton of the outburst of the newly discovered magnetar Swift J1822. 3-1606 (SGR 1822-1606), from 
the first observations soon after the detection of the short X-ray bursts which led to its discovery, 
through the first stages of its outburst decay (covering the time-span from July 2011, until end of 
April 2012). We also report on archival i? 05*^47 observations which witnessed the source during its 
likely quiescent state, and on upper limits on Swift J1822.3-1606's radio-pulsed and optical emission 
during outburst, with the Green Bank Telescope (GBT) and the Gran Telescopio Canarias (GTC), 
respectively. Our X-ray timing analysis finds the source rotating with a period of P = 8.43772016(2) s 
and a period derivative P = 8.3(2) x 10~^^ ss~^ , which entails an inferred dipolar surface magnetic 
field of i? ~ 2.7 x 10^'^ G at the equator. This measurement makes Swift J1822. 3-1606 the second 
lowest magnetic field magnetar (after SGR 0418-1-5729; Rea et al. 2010). Following the flux and 
spectral evolution from the beginning of the outburst, we find that the flux decreased by about an 
order of magnitude, with a subtle softening of the spectrum, both typical of the outburst decay 
of magnetars. By modeling the secular thermal evolution of Swift J1822. 3-1606, we find that the 
observed timing properties of the source, as well as its quiescent X-ray luminosity, can be reproduced 
if it was born with a poloidal and crustal toroidal fields of Bp ~ 1.5 x 10^** G and Btor ~ 7 x 10^** G, 
respectively, and if its current age is ~550kyr. 
Subject headings: stars: magnetic fields — stars: neutron — X-rays: Swift J1822. 3-1606 



1. INTRODUCTION 

The availability of sensitive, large field-of-view X-ray 
monitors such as the Burst Alert Telescope on board 
Swift, and the Gamma-ray Burst Monitor on Fermi, 
makes us witness a golden age for magnetar studies. 
Since the discover y of the first magneta r outbursts 
(IGavriil et al.ll2002l:' Kouveliotou et aLll2003l : iKaspi et al.l 



120031 : llbrahim et al.ll2004l ) . five new members of the class 
have been discovered through the serendipitous detection 
of the typical short X-ray bursts emitted by these highly 
energetic X-ray pulsars, and t he accompanying incre ase 
in the persistent emission (see lRea fc Espositcill20lil for 
a recent review on magnetar outbursts). 

Magnetars, usually recognized in the anomalous X- 
ray pulsar (AXP) and soft gamma-ray repeater (SGR) 



classes, are isolated neutron stars with bright persis- 
tent X-ray emission (Lx ~ lO'^^-lO'^^ ergs~^), rotat- 
ing at spin periods of ^0.3-12 s an d with large pe- 
riod derivatives (10~" -10"^° ss'^ see lMereghettill2008l : 
iRea fc Espositol 120111 for a review) . Sporadically, they 
emit bursts and flares which can last from a fraction of 
seconds to minutes, releasing ^^10^^-10'*^ ergs~^, and are 
often accompanied by long-lived (up to years) increases 
of the persistent X-ray luminosity (outbursts). 

The broadband emission of these objects and their 
flaring activity are believed to be connected to their 
high dipolar and/or toroidal magnetic field: this is in- 
deed supported by the measurement of surface dipo- 
lar B field usually of the order of lO^-lO^^ G (in- 
ferred through the assumption that, as ordinary pul- 
sars, they are spun down via magnetic dipolar losses: 
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Fig. 1. — Left panel: Flux decay of Swift J1822. 3-1606 in the 1-10 kcV energy range (see also TablelSjl. Right panel: pulse phase evolution 
as a function of time, together with the time residuals (lower panel) after having corrected for the linear component (correction to the 
P value). The solid lines in the two panels mark the inferred P-P coherent solution based on the whole dataset, while the dotted lines 
represent the P-P coherent solution based on the data collected during the first 90 days only (see text for the details). 



B - 3.2 X 101^(PP)1/2G, where P is the spin pe- 
riod in seconds, P its first derivative and we assumed 
a neutron star mass and radius of i? ^ 10® cm and 
M ~ 1.4 Mq, respectively). However, the recent detec- 
tion of an SGR showing all the typical emission prop- 
erties defining a ni agnetar (|van der Horst et all 120101 : 
lEsposito et al.l [20To() . but with an inferred dipolar sur- 
face magnetic field < 7.5 x lO^^ Q (jRea et al.ll20Tol ). has 
put into question the need of a high dipolar magnetic 
field (namely higher than the quantum electron critical 
field Bq = meC^/eh ~ 4.4 x 10^^ G ) for an object to 
show magnetar-like activity. 

On 2011 July 14, a new SGR-like burst and 
associated outburst were discovered by the Swift 
Burst Alert Telescope (BAT) , and followed s o on af - 
ter by all X-ray satellites (jCummings et all 1201 1[ ). 
The fast slew of the Swift X-ray telescope (XRT) 
promptly detected a new bright X-ray source at RA : 
18h22'"18!00, Dec : -16°0 4'26"8 (J2000; 1"8 error 
at a 90% confidence level; iPagani et al.l 120 11[). with 
a sp in period of P ~ 8.43 s ()Gogii^ fc Kouveliot"oul 
120 111) . The lack of an optica l /infrared c ounter- 
part (Bandvo padhvavet al.l 120 11: Rc a et al.l l2011bl: 
Ide Ugartc Postigo fc Munoz-Daria&.2011.) . as well as the 
characteristics of the bursts, the X-ray spin period and 
its spe ctral properties (jEsposito et al.l[2011al |b HR"ea et al.l 
l2011al) , led to its identification as a new ma gnetar can- 
didate (jCummings et al.ll201lHHalpernll2(5TTl ). 

After its discovery, many attempts were made to 
measur e the spin period d erivative of Swift J 1 822.3- 
1606 (IGogus et all l20Tlt IKuiper fc HermsenI I20TTI : 
[Livingstone et al.ll2011| ) in order to estimate its surface 



dipolar field. We present here the timing and spectral 
results of the first 9 months of X-ray monitoring (^and 
^ of the new magnetar candidate Swift J1822. 3-1606, 
a detection of its quiescent counterpart in archival data 
(§??), as well as upper limits on its emission in the opti- 
cal and radio bands (^and ^. A detailed study of the 
SGR-like bursts, precise X-ray position and pulse profile 
modeling will be reported elsewhere (Kouveliotou et al. 
in preparation). Using our timing and spectral results, 
we model the source outburst decay, and secular evolu- 
tion, which resulted in an estimate of the its real age and 
crustal toroidal field (§7). 

2. X-RAY OBSERVATIONS AND DATA REDUCTION 

In this study, we used data obtained from different 
satellites (see Table [T] for a summary). Observations and 
data analysis are briefiy described in the following. 

2.1. Swift data 

The X-Ray Telescope (XRT; iBurrows et al|[2005h on- 
board Swift uses a front-illuminated CCD detector sen- 
sitive to photons between 0.2 and 10 keV. Two main 
readout modes are available: photon counting (PC) and 
windowed timing (WT). PC mode provides two dimen- 
sional imaging information and a 2.5073 s time resolu- 
tion; in WT mode only one-dimensional imaging is pre- 
served, achieving a time resolution of 1.766 nis. 

The XRT data were uniformly processed with xrt- 
PIPELINE (version 12, in the heasoft software package 
version 6.11), filtered and screened with standard cri- 
teria, correcting for effective area, dead columns, etc. 
The source counts were extracted within a 20-pixel ra- 
dius (one XRT pixel corresponds to about 2'.'36). For the 



A NEW LOW-B MAGNETAR: Swift J1822.3-1606 



TABLE 1 
Swift J1822. 3-1606 observations used for this work. 



Instrument 


Obs.ID 


Date" 


Exposure 






(MJD TBD) 


(ks) 


Swift/XRT 


00032033001 (PC) 


55 757.75058 


1.6 


RXTE/PCA 


96048-02-01-00 


55 758.48165 


6.5 


Swift/XRT 


00032033002 (WT) 


55 758.68430 


2.0 


Swift/XRT 


00032033003 (WT) 


55 759.69082 


2.0 


RXTE/PCA 


96048-02-01-05 


55 760.80853 


1.7 


Swift/XRT 


00032033005 (WT) 


55 761.54065 


0.5 


RXTE/PCA 


96048-02-01-01 


55 761.55969 


5.0 


Swift/XRT 


00032033006 (WT) 


55 762.24089 


1.8 


RXTE/PCA 


96048-02-01-02 


55 762.47384 


4.9 


Swift/XRT 


00032033007 (WT) 


55 763.30400 


1.6 


RXTE/PCA 


96048-02-02-00 


55 764.61846 


6.1 


RXTE/PCA 


96048-02-02-01 


55 765.46687 


6.8 


Swift/XRT 


00032033008 (WT) 


55 765.85252 


2.2 


Swift/XRT 


00032033009 (WT) 


55 766.28340 


1.7 


RXTE/PCA 


96048-02-02-02 


55 767.59064 


3.0 


RXTE/PCA 


96048-02-02-03 


55 769.35052 


3.4 


Swift/XRT 


00032033010 (WT) 


55 769.49531 


2.1 


Swift/XRT 


00032033011 (WT) 


55 770.39936 


2.1 


Chandra/URC-l 


13511 


55 770.83049 


11.7 


Swift/XRT 


00032033012 (WT) 


55 771.23302 


2.1 


RXTE/PCA 


96048-02-03-00 


55 771.34185 


6.8 


Swift/XRT 


00032033013 (WT) 


55 772.40044 


2.1 


RXTE/PCA 


96048-02-03-01 


55 774.34999 


6.9 


RXTE/PCA 


96048-02-03-02 


55 777.85040 


1.9 


Swift/XRT 


00032051001 (WT) 


55 778.10744 


1.7 


Swift/XRT 


00032051002 (WT) 


55 779.18571 


1.7 


RXTE/PCA 


96048-02-04-00 


55 780.85040 


6.7 


Swift/XRT 


00032051003 (WT) 


55 780.49505 


2.3 


Swift/XRT 


00032051004 (WT) 


55 781.49878 


2.3 


RXTE/PCA 


96048-02-04-01 


55 782.57749 


6.2 


RXTE/PCA 


96048-02-04-02 


55 784.97179 


6.2 


Swift/XRT 


00032051005 (WT) 


55 786.42055 


2.2 


Swift/XRT 


00032051006 (WT) 


55 787.58688 


2.2 


RXTE/PCA 


96048-02-05-00 


55 788.05419 


6.0 


Swift/XRT 


00032051007 (WT) 


55 788.25617 


2.3 


Swift/XRT 


00032051008 (WT) 


55 789.66173 


1.7 


RXTE/PCA 


96048-02-05-01 


55 789.95880 


6.0 


Swift/XRT 


00032051009 (WT) 


55 790.36270 


2.2 


RXTE/PCA 


96048-02-06-00 


55 794.45899 


6.5 


RXTE/PCA 


96048-02-07-00 


55 799.61550 


6.9 


Swift/XRT 


00032033015 (WT) 


55 800.86278 


2.9 


Swift/XRT 


00032033016 (WT) 


55 807.48660 


2.4 


RXTE/PCA 


96048-02-08-00 


55 810.37979 


6.0 


Suzaku/XlS 


906002010 


55 817.92550 


33.5 


RXTE/PCA 


96048-02-10-00 


55 820.23970 


6.7 


Swift/XRT 


00032033017 (WT) 


55 822.82836 


4.9 


Swift/XRT 


00032033018 (WT) 


55 824.71484 


1.5 


RXTE/PCA 


96048-02-10-01 


55826.18540 


5.6 


XMM-Newton 


0672281801 


55 827.25350 


10.6 


Swift/XRT 


00032033019 (WT) 


55 829.45421 


2.3 


Swift/XRT 


00032033020 (WT) 


55 835.54036 


2.6 


RXTE/PCA 


96048-02-11-00 


55835.90370 


7.0 


Swift/XRT 


00032033021 (WT) 


55 842.06040 


4.2 


RXTE/PCA 


96048-02-12-00 


55842.23269 


5.8 


XMM-Newton 


0672282701 


55 847.06380 


25.8 


Swift/XRT 


00032033022 (WT) 


55 849.61916 


3.4 


RXTE/PCA 


96048-02-13-00 


55849.6597976 


5.6 


Swift/XRT 


00032033024 (WT) 


55 862.59155 


10.2 


RXTE/PCA 


96048-02-14-00 


55863.11100 


5.6 


Swift/XRT 


00032033025 (PC) 


55 977.16600 


6.3 


Swift/XRT 


00032033026 (WT) 


55 978.53399 


10.2 


Swift/XRT 


00032033027 (PC) 


55 981.99499 


11.0 


Swift/XRT 


00032033028 (WT) 


55 982.96299 


7.0 


Swift/XRT 


00032033029 (WT) 


55 985.17799 


7.0 


Swift/XRT 


00032033030 (WT) 


55 985.55000 


7.0 


Swift/XRT 


00032033031 (WT) 


55 991.09231 


6.7 


XMM-Newton 


0672282901 


56022.95692 


26.9 


Swift/XRT 


00032033032 (WT) 


56 031.141159 


4.3 



spectroscopy, we used the spectral redistribution matri- 
ces in CALDB (20091130; matrices version vOlS and v014 
for the PC and WT data, respectively), while the ancil- 
lary response files were generated with xrtmkarf, and 
they account for different extraction regions, vignetting 
and point-spread function corrections. 

2.2. RXTE data 



Th e Proportional Counter Array (PCA; iJahoda et alj 
|1996[ ) on-board RXTE consists of five coUimated 
xenon/methane multianode Proportional Counter Units 
(PCUs) operating in the 2-60 keV energy range. Raw 
data were reduced using the ftools package (version 
6.11). To study the timing properties of Swift J1822.3- 
1606, we restricted our analysis to the data in Good 
Xenon mode, with a time resolution of 1 /.ts and 256 en- 
ergy bins. The event-mode data were extracted in the 2- 
lOkeV energy range from all active detectors (in a given 
observation) and all layers, and binned into light curves 
of 0.1 s resolution. 

2.3. Suzaku data 



Suzaku (jMitsuda et all I2007D observed the field of 
Swift J1822.3-1606 on 2012 September 13-14 with the 
pulsar located at the X-ray Imaging Spectrometer (XIS; 
iKovama et al.lf20b7i ) nominal position. The XIS consists 



of three front-illuminated (FI) CCD cameras (XISO, XIS2 
and XIS3), and one that is back-illuminated (BI; XISl). 
One of the FI CCDs, XIS2, was not available at the time 
of our observation. XISl and XIS3 were operating in 
Normal Mode without any option (all the pixels on the 
CCD are read out every 8 s), while XISO was operating 
with the 1/8 Window option allowing a read out time of 
Is. 

For each XIS, 3x3 and 5x5 edit modes cleaned event 
data were combined. Following standard practices, we 
excluded times within 436 s of Suzaku passing through 
the South Atlantic Anomaly and we also excluded the 
data when the line of sight was elevated above the Earth 
limb by less than 5°, or less than 20° from the bright- 
Earth terminator. Moreover, we excluded time windows 
during which the spacecraft was passing through a cut- 
off rigidity of below 6 GV. Finally, we removed hot and 
fiickering pixels. The resulting total effective exposure 
was '^33.5 ks for each XIS. The SGR net count rates are 
0.710(5), 1.180(6), and 1.060(6) count s"! in the XISO, 
XISl and XIS3, respectively. For the spectral analysis, 
we used only XISO and XIS3, which are the best cali- 
brated cameras, while for the timing analysis we made 
use only of the XISO data, which owing to the 1/8 Win- 
dow option, have a timing resolution adequate to sample 
the pulsar spin period. 

2.4. Chandra data 

The Chandra X-ray Observatory has observed 
Swift J1822. 3-1606 with the High Resolution Imaging 
Camera (HRC~I; Zombeck et al. 1995) on 2011 July 
28, for - 12ks (ObsID: 13511). Data were analyzed 
using standard cleaning proceduretQ and CIAO version 
4.4. Photons were extracted from a circular region with 
a radius of 3" around the source position, including 

^ http://asc.harvard.edu/ciao/threads/index.html 



Mid-point of the observations. 
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Fig. 2. — Spectra of the Suzaku (blue) and XMM-Newton observations fitted together with a blackbody plus power-law and two-blackbody 
models (see text and Table[2]for details). 



TABLE 2 
Spectral analysis of the Suzaku and XMM-Newton data. 









Blackbody + Powerlaw 




Two Blackbodies 




Instrument 


Time" 


Flux'' 


kT (keV) 


Rbb (km)= r 


kTi (kcV) 


^BBi (km)= kTa (keV) 


flBB2 (km)'= 


Suzaku 


60.93±0.48 


1.78±0.01 


0.678±0.008 


I.2±0.f 2.90±0.04 


0.39±0.0I 


2.6±0.4 0.79±0.0f 


0.9±0.1 


XMM-Newton 


70.25±0.06 


1.70±0.01 


0.689±0.006 


f.liO.f 2.86±0.03 


0.40±0.0I 


2.6±0.4 0.84±0.0f 


O.SitO.l 


XMM-Newton 


90.06±0.11 


1.20±0.03 


0.679±0.005 


l.OiO.l 2.99±0.03 


0.37±0.02 


2.5±0.4 0.79±0.0f 


0.7±0.1 


XMM-Newton 


266.f±0.2 


0.4±0.1 


0.623±0.008 


0.6f±0.08 3.05±0.04 


0.35±0.01 


2.0±0.2 0.78±0.0f 


0.4±0.1 



°' Times arc calculated in days from MJD 55 757.0 . 

Fluxes are in units of 10~ ergs" cm~ , referred to the BB+PL fit, and calculated in the 1— lOkeV energy range. Errors in the table are 
given at la confidence level. Reduced x^ and absorption values are x^/dof =1.05/2522 and A^h = 0.50(1) X lO^^cm"^, x^/dof =1.06/2522 
and A^H = 0.21(1) X lO^^cm"^, for the BB+PL and BB+BB models, respectively. 



^ Radii are calculated assuming a distance of 5kpc. 

more than 90% of the source photons (see Kouvehotou 
et al. 2012 in prep for further details on this obser- 
vation). We inferred an effective HRC-I count-rate of 
70 ±1 counts s"^ 

2.5. XMM-Newton data 

We observed Swift J1822 3-160 6 three time with 
XMM-Newton pansen et all |2001[ ) on 2011 September 
23, October 12, and April 05, for ~10, 25, and 27 ks, 
respectively. Only the second observation was partially 
affected by background flares which we have cleaned dur- 
ing our spectral analysis (see §3.1|) resulting in a net ex- 
posure time of 19.4 ks. Data have been processed us- 
ing SAS version 11, and we have employed the most up- 
dated calibration files available at the time the reduction 
was performed (April 2012). Standard data screening 
criteria are applied in the extraction of scientific prod- 
ucts. For our analysis we used only the EPIC-pn camera, 
and we checked that the two MOS cameras gave con- 
sistent results. The EPIC-pn camera was set in Prime 
Large Window mode (timing resolution 48 ms) , with the 
source at the aim-point of the camera. We have extracted 



the source photons from a circular region of 30"radius, 
and a similar region was chosen for the background in 
the same CCD but as far as possible from the source 
position. We restricted our analysis to photons hav- 
ing PATTERN<4 and FLAG=0. We find a EPIC-pn 
(background-subtracted) count rate of 5.03(3), 3.68(2) 
and 1.42(1) counts s~^ for the first, second, and third ob- 
servation, respectively. 

3. RESULTS OF THE X-RAY MONITORING 

3.1. X-ray spectral modeling 

Spectra were extracted as explained in ^ and re- 
binned in order to have at least 20 counts per bin in 
the Swijt spectra, and 50 counts per bin in the Suzaku 
and XMM-Newton spectra. We started our spectral 
analysis by fitting our higher quality spectra, those 
from the three XMM-Newton/ pn, and Suzaku/XISOS ob- 
servations, with several models (using XSPEC version 
12.7.0; see Figured and Table [2]). Best fits were found 
using a blackbody plus power-law (BB-t-PL; xj/dof 
= 1.05/2522) and a 2 blackbodies (2BBs; xl/doi = 
1.06/2522) model, all corrected for the photoelectric ab- 
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Fig. 3. — Left panel: Spectral parameters for a blackbody plus power-law fit. Right panel: Spectral parameters for a two-blackbody fit 
(see text and Table[3]for details). 



sorption (phabs model with solar abundances assumed 
from Anders & Grevesse (1989) and photoelectric cross- 
section from Balucinska-Church & McCammon (1998)). 
The hydrogen column density along the line of sight 
was fixed to the same value for all of the spectra for 
a given model. We obtained Nh = 0.50(1) and 0.21(1) x 
lO^^cm"^ (errors in the text are at la level unless other- 
wise specified), for the BB+PL and 2BBs model, respec- 
tively (see Table [2|). In Figure [2] we show the residuals of 
this spectral modeling, and note that, although statisti- 
cally the fits can be considered equally good, the BB+BB 
model departs from the data at higher energies. 

Already from this first analysis it is evident how the 
spectrum is changing in time, although very slowly. 

We then expanded our spectral modeling by fitting 
simultaneously all the Swift/XRT, S'wzafcu/XISOS, and 
XMM-Newton/pn spectra. Again both models gave sat- 
isfactory fits (see Table[3]). The hydrogen column density 
along the line of sight was fixed to the same values found 
from the modeling of our previous analysis. Table|3]sum- 
marizes the obtained spectral parameters. 

For the BB-I-PL model, we first allowed all parameters 
to vary freely, and we noticed that the BB temperature 
was consistent with being constant in time in the early 
phases of the outburst (most probably changing too little 
for our spectral analysis to be sensible to its variations). 
This was visible already from Table [2] when consider- 
ing only the most detailed spectra. We then tied the BB 
temperature across all spectra in the first 100 days of the 
outburst, and similarly we did for the last spectra (be- 
tween 200-300 days after the trigger). Best fit (reduced 
X^/dof = 1.1/6501) was found with a BB temperature 
of kT = 0.660(8) keV and kT = 0.621(7) keV, for the 
early and late times spectra, respectively. More detailed 
spectra would have certainly disentangled a slow decay 
between those two values. Figure [3] (left panel) shows 
the time evolution of the power-law index (F) and the 
BB area. We can see how the latter shrinks as the out- 
burst decays, while the power law index increases slowly, 
anti-correlated with the X-ray flux (see also FigureH] left 
panel) . 



On the other hand, for the BB+BB model we noticed 
that, by leaving all the parameters free to vary across 
all the spectra, the temperature and the radius of the 
first blackbody were not varying significantly in time. 
Similarly to the BB-I-PL case, we then fixed those values 
to be the same in all spectra at early and late times 
separately. This resulted in the best fit values of kTi = 
0.388(8) keV and BB Radiusi = 2.5(1) km, and kTi = 
0.358(7) keV and BB Radiusi = 1.9(1) km (to estimate 
the BB radii we assume a source distance of 5kpc). The 
best fit had a reduced x^/dof = 1.1/6542. The second 
blackbody has a relatively steady temperature around 
^0.7 keV (see Fig. [3] right panel) and its radius shrinks 
during the outburst decay. 

In the late time spectra, taken >200 days after the 
outburst onset, the source fiux decreases substantially 
(from ~ 24 to 0.4 x lO^^^ergs^^cm^^), as the spectrum 
continues to soften. 

We have also tried to model t he spectra with a reso- 
nant cyclotron s cattering model (jRea et al.l 120071 120081 : 
iZane et al]l2009l ). and, although the fits gave a good chi- 
squared value {xt '^l. 1/5 912), the low magnetic field 
of Swift J1822.3-1606 (see ^321) makes the use of those 
models, envisaged for B ^ 10^** G, questionable (see also 
TuroUa et al. 2011). 

We note that, although the x^ values of both the 
BB+PL and BB+BB fits might appear not acceptable 
from a purely statistical point of view, many systematic 
errors are present in the simultaneous spectral modeling 
of different satellites (the most severe being the uncer- 
tainties in the inter-calibration between them, which is 
believed to be within a 5% error). We did not add any 
systematic error in the spectral fitting to show the pure 
residuals of the fit; however, with only 5% systematic er- 
ror, the reduced x^ values would decrease substantially, 
reaching a fully acceptable fit for both models {xt '^1.0). 

3.2. X-ray timing analysis 

For the X-ray timing analysis we used all data listed 
in Table[l] after referring the event arrival times to the 
barycenter of the Solar System (assuming the source co- 
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TABLE 3 

Spectral analysis 









Blackbody + Powerlaw* 


Two Blackbodies" 


ObsID 


Time"^ 


Flux*" 


r 


BE radius {km)'= 


kT2 (keV) 


BB2 radius (km)'= 


33001 


0.76±0.01 


23.4±1.0 


0.4±1.0 


4.8±0.4 


0.S0it0.02 


3.37it0.13 


33002 


1.69±0.01 


24.1±0.5 


2.0±0.1 


5.4±0.1 


0.72it0.01 


4.7it0.1 


33003 


2.70±0.01 


18.8±0.6 


2.3±0.1 


4.6±0.1 


0.70it0.01 


4.4it0.1 


33005 


4.11±0.01 


18.8±0.6 


2.4±0.1 


3.3±0.4 


0.68it0.02 


4.4it0.2 


33006 


5.25±0.20 


15.8±0.5 


2.4±0.1 


4.0±0.1 


0.69it0.01 


4.2it0.1 


33007 


6.31±0.08 


17.8±0.7 


2.1±0.1 


4.2±0.1 


0.71it0.01 


4.1it0.1 


33008 


8.86±0.07 


12.5±0.5 


2.2±0.1 


3.8±0.1 


0.71it0.01 


3.5it0.1 


33009 


9.29±0.04 


12.2±0.4 


2.4±0.1 


3.4±0.1 


0.69it0.01 


3.6it0.1 


33010 


12.50±0.04 


11.2±0.4 


2.1±0.1 


3.3±0.1 


0.73it0.01 


3.1it0.1 


33011 


13.40±0.07 


10.2±0.3 


2.4±0.1 


3.4±0.1 


0.70it0.01 


3.2it0.1 


33012 


14.24±0.04 


10.0±0.5 


2.2±0.2 


3.1±0.1 


0.71it0.01 


3.0it0.1 


33013 


15.27±0.01 


9.3±0.1 


2.2±0.1 


2.3±0.1 


0.71it0.01 


2.8it0.1 


51001 


21.11±0.10 


7.0±0.4 


2.2±0.1 


2.2±0.1 


0.74it0.01 


2.3it0.1 


51002 


22.19±0.16 


6.3±0.3 


2.4±0.1 


2.0±0.2 


0.70it0.01 


2.5it0.1 


51003 


23.50±0.08 


6.3±0.2 


2.4±0.1 


2.0±0.1 


0.70it0.01 


2.4it0.1 


51004 


24.50±0.14 


6.1±0.3 


2.3±0.1 


1.5±0.2 


0.68it0.01 


2.4it0.1 


51006 


30.59±0.01 


4.5±0.2 


2.6±0.2 


1.9±0.2 


0.68it0.01 


2.2it0.1 


51007 


31.07±0.07 


4.1±0.2 


2.7±0.2 


2.0±0.1 


0.69it0.02 


2.1it0.1 


51008 


31.40±0.01 


4.0±0.2 


2.5±0.1 


1.9±0.1 


0.72it0.01 


l.SitO.l 


51009 


32.67±0.01 


3.5±0.3 


2.9±0.1 


1.8±0.2 


0.67it0.01 


2.0it0.1 


33015 


43.87±0.10 


2.7±0.3 


2.8±0.2 


l.SitO.l 


0.71it0.02 


l.SitO.l 


33016 


50.49±0.31 


2.4±0.4 


3.0±0.3 


1.6±0.1 


0.71it0.02 


1.4it0.1 


6002010 


60.93±0.48 


1.78±0.03 


2.85±0.03 


1.30±0.01 


0.77it0.01 


1.03it0.02 


33017 


65.83±0.13 


1.7±0.2 


2.6±0.3 


1.2±0.1 


0.78it0.02 


0.94it0.05 


33018 


67.72±0.20 


2.1±0.2 


2.3±0.2 


l.OitO.l 


O.Slit0.04 


l.OitO.l 


81801 


70.25±0.06 


1.70±0.03 


2.76±0.03 


1.19±0.01 


0.78it0.01 


0.96it0.02 


33019 


72.45±0.35 


l.SitO.l 


2.3±0.6 


l.SitO.l 


0.70it0.02 


1.3it0.1 


33020 


78.54±0.46 


l.SitO.l 


2.7±0.2 


l.OitO.l 


0.71it0.03 


l.litO.l 


33021 


85.06±0.38 


1.48±0.15 


2.4±0.6 


1.4it0.1 


0.77it0.02 


0.94it0.06 


82701 


90.06±0.11 


1.20±0.03 


2.91±0.03 


1.03it0.01 


O.S3it0.01 


0.67it0.02 


33022 


92.62±0.38 


1.4±0.1 


3.1±0.2 


1.2it0.1 


0.75it0.02 


0.92it0.07 


33024 


99.57±0.37 


0.98±0.06 


2.8±0.2 


O.SitO.l 


0.70it0.02 


O.S7it0.04 


33025 


220.2±0.8 


0.55±0.08 


2.6±0.2 


O.SitO.l 


O.SOitO.05 


0.5it0.2 


33026 


221.5±0.4 


0.54±0.06 


2.9±0.3 


O.SitO.l 


0.74it0.04 


0.6it0.2 


33027 


225.0±0.1 


0.47±0.08 


2.9±0.4 


0.9it0.1 


0.78it0.05 


0.5it0.2 


33028 


226.0±0.1 


0.47±0.08 


3.4±0.5 


O.SitO.l 


0.72it0.05 


0.6it0.2 


33029 


228.2±0.1 


0.46±0.08 


3.0±0.5 


0.9it0.1 


0.73it0.05 


0.6it0.2 


33030 


228.5±0.2 


0.54±0.08 


2.5±0.3 


o.eio.i 


O.SOitO.05 


0.5it0.2 


33031 


234.5±0.4 


0.50±0.08 


3.1±0.5 


O.SitO.l 


0.72it0.08 


0.6it0.2 


82901 


266.1±0.2 


0.41±0.02 


3.06±0.04 


0.63it0.02 


0.78it0.01 


0.45it0.08 


33032 


274.5±0.4 


O.SitO.l 


2.8±0.4 


o.eio.i 


0.75it0.06 


0.6it0.2 



* The absorption value and the blaekbody temperature were fixed to be the same for all spectra of the first and second set: A^h — 0.50(1) X 
lO^^cm"^, kT = 0.660(8) keV and kT = 0.621(7) keV, for the first and second set of spectra, respectively. Reduced xj/dof = 1.1/6501. 

*'* The absorption value, blackbody temperature and radius were fixed to be the same for all spectra of the first and second set: A^h — 
0.21(1) X lO^^cm-^, fcTi = 0.388(8) keV and BBi Radius = 2.5(1) km, and kTi = 0.358(7) keV and BBi Radius = 1.9(1) km, for the first 
and second set of spectra, respectively. Reduced x^/dof — 1.1/6542 . 

^ Times are calculated in days from MJD 55 757.0 . 

Fluxes are in units of 10^ ergs^ cm^ , referred to the BB+PL fit, and calculated in the 1— lOkeV energy range. Errors in the table are 
given at Icr confidence level. 

^ Radii are calculated assuming a distance of 5kpc. 



ordinates bv lPagani et al.ll2011l and the DE200 ephemeris 
for the Solar System). The first Swift/ XRT event Hsts 
were used in order to start building up a phase coherent 
timing solution and to infer the SGR timing properties. 
We started by obtaining an accurate period measurement 
by folding the data from the first two XRT pointings 
which were separated by less than Iday, and studying 
the phase evolution within thes e observations by me ans 
of a phase- fitting technique fsee lDaU'Osso et "aLll2003l for 



details). Due to the possible time variability of the pulse 
shape we decided not to use a pulse template in the cross- 
correlation, which might artificially affect the phase shift, 
and we instead fit each individual folded light curve with 
two sine functions, the fundamental plus the first har- 
monic. In the following we also implicitly assume that 
the pulsation period (and its derivative) is a reliable esti- 
mate of the spin period (and its derivative), an assump- 
tion which is usually considered correct for isolated neu- 
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Fig. 4. — Pulse profiles with 32 bins and relative pulsed fractions as a function of the energy band, for the three XMM-Newton observations. 



tron stars. 

The rcsufting best-fit period (reduced x^ = 1-1 for 2 
dof) is P = 8.43966(2) s (all errors are given at la c.l.) 
at the epoch MJD 55757.0 . The above period accu- 
racy of 20 /is is enough to phase-connect coherently the 
later Swift, RXTE, Chandra, Suzaku, and XMM-Newton 
pointings. The procedure was repeated by adding, each 
time, a further observation folded at the above period, 
and following the phase evolution of the ascending node 
of the fundamental sine function best fitting the pro- 
file of each observation. The relative phases were such 
that the signal phase evolution could be followed unam- 
biguously for the whole visibility window until November 
2011 (see Figure[l]). When adding the RXTE dataset, we 
also corrected the output phases by a small constant off- 
set (^0.02), likely due to the different energy ranges and 
responses. 

We modeled the phase evolution with a polynomial 
function with a linear plus quadratic term, the inclu- 
sion of the latter results in a significant improvement 
of the fit (an F-test gives a probability of 7 x 10^^ 
that the quadratic component inclusion is not required) . 
The corresponding coherent solution (valid until Novem- 
ber 2011) is P = 8.43772007(9) s and period derivative 
P = 1.1(2) X 10-13 s s-i (x^ = 132 for 57 dof; at epoch 
MJD 55757.0). The above solution accuracy allows us 
to unambiguously extrapolate the phase evolution until 
the beginning of the next Swift visibility window which 
started in February 2012. 

The final resulting phase-coherent solution (see also 
TableHJ, once the latest 2012 observations are included, 
returns a best-fit period of P = 8.43772016(2) s and 
period derivative of P = 8.3(2) x 10"" s s^^ at MJD 
55757.0 (y^ = 145 for 67 do f; preliminary results were 
reported in [Israel et al.|[20T^ . The above best-fit values 
imply a surface dipolar magnetic field of P ~ 2.7 x lO^^^ G 
(at the equator), a characteristic age of Tc = P/2P ~ 
1.6 Myr, and a spin-down power Lrot — AnlP/P^ ~ 
1.7x10^" erg s-i (assuming a neutron star radius of 10 km 
and a mass of 1.4Mq). 

The final solution has a relatively high r.m.s. (^ 
120 ms) resulting in a best-fit reduced xt — 2.1. The 
introduction of high-order period derivatives in the fit of 
the phase evolution does not result in a significant im- 
provement of the fit (x^ = 135 for 66 dof; F-test gave 
a probability of 0.03 that the cubic component inclusion 
is not required) . This results in a 3cr upper limit of the 



second derivative of the period of P < 5.8 x 10 ^^s s ^ . 

These values of P and P are in agreement (within la) 
with those inferred for the 2011 visibility window re- 
ported abov43 However, this solution is not consistent, 
within 3a, with those already reported in the litera- 
ture and based on a reduced d ataset (Kui per fc Hermseiil 
120111: [Livingstone et al.|[20llt valid until 70 and 90 days 
from the onset of the outburst, respectively). In order to 
cross check our results and to compare them with those 
previously reported we fit only those observations of our 
dataset within about 90 days from the trigger. The corre- 
sponding best-fit parameters are P = 8.4377199(1) s and 
period derivative of P = 1.6(4) x lO'^^s s"^ (x^ = 119 
for 52 dof; at epoch MJD 55757.0). The latter values are 
consistent with those of Livingstone et al. (2011). 

This analysis together with the relatively high r.m.s. 
value suggest that the timing parameters of the pulsar 
are "noisy". Correspondingly, a timing solution based 
on a longer baseline may decrease the effect of a noisy 
behavior, while those reported earlier are likely affected 
by the shorter time-scale variability of the timing param- 
eters. 

In FigurelH we show the pulse profiles as a function 
of energy for the three XMM-Newton observations 
(see ij2.5p . folded with the best- fit timing solution 
reported above. We derived pulsed fractions (defined as 
PF = [Max(counts/s) — Min(counts/s)]/[Max(counts/s)-|- 
Min(counts/s)]) in the 0.3-1, 1-2, 2-3, 3-5, and 5- 
lOkeV energy bands, of 47 ± 2, 56 ± 1, 65 ± 1, 81 ± 1, 
and 86 ± 2 %, for the first observation; for the second 
observation: 50 ± 2, 56 ± 1, 62 ± 1, 67 ± 2, and 81 ± 2 % 
; and for the third observation: 55 ±3, 62 ± 3, 67 ± 3, 
64 ± 3, and 73 ± 4% . 

4. ROSAT PRE-OUTBURST OBSERVATIONS 

The Rontgensatellit [ROSAT) Position Sensitive Pro - 
portional Counter (PSPC; Pfeffe rmann et al.l Il987| ) 
serendipitously observed the region of the sky including 
the position of Swift J1822. 3-1606 between 1993 Septem- 
ber 12 and 13 (obs. ID: rp500311n00), for an effective 
exposure time of 6.7 ks. 

By means of a sliding cell source-detection algorithm, 
we found a 5.5it significant source with a 0.1-2.4 keV 

■^ Note that the two data sets are not independent. However, 
we checked that when deriving a timing solution independently 
using only RXTE, Chandra, XMM-Newton and Suzaku data for 
the first 120 days, and all the Swift observations plus the latest 
XMM-Newton observation for the whole ~300 days, the two (now 
independent) solutions are still in agreement within la. 
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TABLE 4 
Timing parameters for Swift J1822. 3-1606. 



18"22'"48 



Reference Epoch (MJD) 


55757.0 


Validity period (MJD) 


55757- 56032 


P(s) 


8.43772016(2) 


P(ss-i) 


8.3(2) X 10-1-* 


P(sS-2) 


< 5.8 X 10"^^ 


:/(Hz) 


0.118515426(3) 


z>(s-^) 


1.17(3) X 10"^^ 


z>(Hzs-2) 


< 8.1 X 10"^^ 


xVdof 


145/66 


RMS residuals (ms) 


120 


B (Gauss) 


2.7 X 10" 


Lrot (erg s"^) 


1.7 X 10^° 


Tc (Myr) 


1.6 



count rate of 0.012(3)counts s -'^ at the coordinates 
RA = 18^22™18!1, Decl. = -16°04'26'.'4 (positional un- 
certainty of 30" radius at 90% c.L; J2000). This source 
is also listed in the WGA and RXP catalogs, namely 
IWGA J1822.2-1604 and 2RXP J182217.9-160417, with 
consistent values of count rate. The positions of the lat- 
ter objects are ~20" and 10" from the Swift-XRT po- 
sition. Given the relatively large ROSAT/PSPC posi- 
tional uncertainty, we believe the latter two sources and 
Swift J1822. 3-1606 are the same object, which we pro- 
pose as the SGR quiescent counterpart (see Figure[Q. 

We downloaded the relevant files of the ROSAT 
pointed observation and extracted the photon arrival 
times from a circle of ^ 100" radius (corresponding 
to an encircled energy of >9G%) around the X-ray po- 
sition. We found that an absorbed blackbody with 
7Vh< 7 X 10^1 cm-2 (see also iS3|), kT = 0.20 ± 0.05 
keV, and a radius of 5 ± 2 km, best fit the data (reduced 
X^ = 0.8 for 4 degree of freedom). We infer an observed 
flux of « 1.3 X 10^13 erg cm'^ s^^ and « 4 x 10"" 
erg cm~^ s~^ in the 0.1-2.4 keV and 1-lOkeV energy 
ranges, respectively. Assuming a distance of 5kpc, this 
flux results in a bolometric luminosity during the quies- 
cence state of Lqui ~ 4 x lO'^^ergs"^. 

No significant periodic signal was found by means of 
a Fourier Transform, even restricting the search around 
the 8.44 s period. The 3cr upper limits on the pulsed 
fraction (semi-amplitude of the sinusoid divided by the 
source average count rate) is larger than 100%. 

5. OPTICAL AND INFRARED OBSERVATIONS 

We observed the field of Swift J1822.3-1606 with the 
10.4-m Gran Telescopio Canarias (GranTeCan) at the 
Roque de los Muchachos Observatory (La Palma, Spain). 
Images were taken in service mode on 2011 July 21 with 
the OSIRIS camera, a two-chip Marconi CCD detector 
with a nominal un- vignetted 7.'8 x 7.'8' field of view, and 
an unbinned pixel size of 0"125. Observations were taken 
through the Sloan z filter (A = 969.4 nm; AA — 261 nm). 
We used a 5-point dithering pattern to correct for the ef- 
fects of the CCD fringing in the Red part of the spectrum. 
We accurately selected the pointing of the telescope to 
position our target in the right CCD chip and a bright 
{B ~ 10) star ~54" East of it in the left one, in order to 




16°12 00 



Fig. 5.— ROSAT 1993 image of the region of Swift J1822. 3-1606. 
The circle is centred at the Swift/XHT position. 

avoid the contamination from ghost images and satura- 
tion spikes. Unfortunately, the observations were taken 
in conditions of very high sky background due to the high 
lunar illumination, with the Moon phase at ~0.5 and 
angular distance < 90°, and with a seeing ranging from 
l-2"5. Observations were taken using exposure times of 
108 and 54 s, with the latter chosen to minimize the sky 
background induced by the Moon. The total integration 
time was 4100 s. We reduced the images with the dedi- 
cated tools in the IRAF ccdred package for bias subtrac- 
tion and flat-fielding, using the provided bias and sky flat 
images. We performed the photometric calibration using 
exposures of the standard star PG 15284-0628. In order 
to achieve the highest signal-to-noise, we filtered out 
observations taken with the highest seeing and sky back- 
ground. We aligned and co-ad ded all the best-im ages 
by means of the swarp program (iBertin et al.ll2002l ). ap- 
plying a 3cr filter on the single pixel average to filter out 
residual hot and cold pixels and cosmic ray hits. We per- 
formed the astrometry calibration of the OSIRIS image 
with the WCStools astrometry packagcQ, using as a ref- 
erence the coordinates of sta rs selected from the GSC2 
catalogue (jLasker et al.ll2008l ). Due to the significant and 
unmapped CCD distortions, we only obtained an rms of 
0"3 on the astrometric fit. We detected three objects (SI, 
S2, S3) wi thin or close to th e Swift J1 822.3-1606 po sition 
(see also iRea et all 1201 Ibt fCorosabel et all [201ll ) . We 
computed their flux through standard aperture photom- 
etry using the IRAF package apphot. Their z-band mag- 
nitudes are 18.13±0.16, 20.05±0.04, and 19.94±0.04, for 
SI, S2 and S3, respectively (see Figure[5|). We detected 
no other object consistent wit h the refined Swif t/XKT 
position of Swift J1822. 3-1606 (|Pagam et al.ll2011[ ) down 
to a 3cr limiting magnitude of z = 22.2 ±0.2. Given their 
bright optical magnitudes, we doubt that any of these 
objects is the optical counterpart to Swift J1822. 3-1606. 
Based upon GranTecan spectroscopy, de Ugarte Postigo 
& Munos-Darias (2011) suggest that SI and S2 are G to 
M-type stars. 
As a reference, we inspected images of the 

^ http://tdc-www.harvard.edu/wcstools/ 
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TABLE 5 

GBT RADIO OBSERVATIONS OF SwiFT J1822. 3-1606. 



Date 


Start time (MJD) 


Exposure (s) 


Smin (mjy)* 


2011-07-22 


55764.26856481 


1031.8269 


0.06 


2011-08-18 


55792.00894675 


967.3408 


0.06 


2011-09-20 


55824.04596064 


1365.0146 


0.05 


2011-10-19 


55853.87988425 


1375.7644 


0.05 



Fig. 6.— Left panel: GranTeCan image of the Swif t J1822.3-1606 field (30" X 30") taken with the OSIRIS camera (z band). The circle 
corresponds to the Swift/XICC position of the source l lPagani et al.|[2blll ). The three objects detected close to, or within, the X-ray source 
error circle are labelled. Right panel: UKIRT image of the same area obtained with the WFC in the K band. 

a month one ft-om the other (see Table [5|). Data were ac- 
quired with the Gree n Bank Uhimate Puls ar Processing 
Instrument (GUPPI; iDuPlain et"all I2008D at a central 
frequency of 2.0 GHz over a total observing bandwidth 
of 800 MHz. 

For each observation, in order to correct for the dis- 
persive effects of the interstellar medium, the bandwidth 
was split into 1024 channels about 250 of which were un- 
usable because of radio frequency interferences (RFI), 
leaving us with 600 MHz of clean band. The obser- 
vations lasted 16 to 23 minutes and were sampled ev- 
ery 0.6557 ms. Since the pulsar rotational parameters 
are known from X-ray observations (see i j3.2p . we first 
folded the data at the known period. We also folded 
the data at half, one third and a quarter of the nom- 
inal period in order to detect putative higher harmon- 
ics components of the intrinsic signal, in case the latter 
were deeply contaminated by RFI. Folding was done us- 
ing DSPSR (van Straten & Bailes 2011) to form 30 s long 
sub-integrations subdivided into 512 time bins. The sub- 
integrations and the 1024 frequency channels, cleaned 
from RFI, were then searched around the nominal period 
P and over a wide range of dispersion measure (DM) val- 
ues (from to 1000 pccm^'^) to find the P-DM combina- 
tion maximizing the signal-to-noise ratio. No dispersed 
signal was found in the data down to a flux of about 
0.05 mJy depending on the observation (see Tabled. 

Data were also blindly searched using the code suites 
PRESTC0 and SIGPROCO- In both cases, after de- 
dispersion of the data with 839 trial DMs (ranging from 
to 1000 pccm^'^) and removal of the frequency chan- 
nels affected by RFI, the time series are transformed 
using fast Fourier algorithms and their power spectra 
searched for relevant peaks. These Fourier-based search 
techniques require a 2" number of time samples in input; 
for this reason the amount of data analyzed was 1030 s, 
(a minute of fake data were added to the shortest ob- 
servation) about two-thirds of the total of the longest 
observation, hence the flux limit attained, depending on 
the inverse of the square root of the integration time, was 
proportionally higher. With SiGPROC we also searched 



* S„iin is the minimum flux density reached. 

Swift J1822. 3-1606 field taken prior to our GranTeCan 
observations, i.e. when Swift J1822. 3-1606 was probably 
in quiescence. To this aim, we searched for near-infrared 
(IR) observations taken as part of the UK Infrared Deep 
Sky Survey (UKIDSS; Lawrence et al. 2Q 0l), performed 
with the Wide Field Camera (WFCAM; (Casali et al] 
120071 ) at the UK Infrared Telescope (UKIRT) at the 
Mauna Kea Observatory (Hawaii). The Swift J1822.3- 
1606 field is indeed included in the UKIDSS Galactic 
Plane Survey (GPS) and data are available through Data 
Rele ase 8 plus. Observations we re taken on 2006 May 
3rd (jBandvopadhvav et al.ll20l"TI ). We downloaded the 
fully reduced, calibrated, and co-added J, H, i^-band 
science images of the Swi ft J1822. 3-1606 field produced 
by the UKIDSS pipeline (iHamblv et al.ll2008[ ) together 
with the associated object catalogues through the WF- 
CAM Science Archive (WSA0 in terface. The WFCAM 
astrometry is based on 2MASS (iSkrutskie et alJl2006D 
and is usually accurate to ^ O'.'l (jLawrence et al.ll2007[ ). 
We clearly identified objects SI (J = 13.92; H = 12.37; 
K = 11.62), S2 (J = 16.62; H = 15.75; K = 15.20), S3 
(J = 16.43; H ^ 15.40; K = 14.88) in the UKIDSS im- 
ages (see Figure[5|) , with a relative fiux comparable to the 
z-band flux measured on the OSIRIS ones. No other ob- 
ject is detected at the Swift J1822. 3-1606 position down 
to 5a limiting magnitudes of J ^ 19.3, H '^ 18.3 and 
K ^ 17.3. 

6. RADIO OBSERVATIONS 

Radio observations of Swift J1822. 3-1606 were per- 
formed at the 101 m Green Bank Telescope (GBT) on 
four occasions after the X-ray outburst, spaced by about 

* http://survcys.roe.ac.uk/wsa/ 



^ See http : //www ■ cv . nrao . edu/$\siin$sr ansom/presto/ [ 
^ See http://sigproc.sourccforge.net/ 
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Fig. 7. — From top left to bottom right, the evolution of the period, period derivative, luminosity, and surface dipole field, according to 
the model discussed in the text. Grey lines report on the current values of the Swift J1822. 3-1606 period and period derivative (with its 
error). The vertical dashed lines mark the source "real" age of 550 kyr . 



the data for single de-dispersed pulses but no signal was 
found in either the Fourier domain or the single pulse 
searches. 

7. DISCUSSION 

7.1. The secular thermal evolution of 
Swift J 1822.3-1606 

To investigate whether the observed properties of 
Swift J1822. 3-1606 are consistent with those of an 
evolved magnetar, as suggested by its characteristic age 
of Tc ~ 1.6 Myr, we followed the secular evolution of this 
object using a two-dimensional magneto-thermal evolu- 
tion code. We refer to Pons et al. (2009) and Aguilera 
et al. (2008) for details about the code and the micro- 
physical inputs of the model. This allows us to esti- 
mate the natal properties of the neutron star, its cur- 
rent age and internal field strength. We have considered 
the evolution, including magnetic field decay and heating 
by Ohmic diffusion, of an ultra-magnetized neutron star 
with a mass of 1.4 M0, with no exotic phases nor fast 
neutrino cooling processes, but with enhanced neutrino 
emission from the breaking and formation of neutron and 
proton Cooper pairs (standard cooling scenario). We as- 
sumed an initial neutron star spin period of 10 ms and an 



initial dipolar field of Bdip(i = 0) = 1.5 x lO^"* G. In Fig- 
ure [7] we plot the evolution of spin period, period deriva- 
tive, luminosity, and the dipolar surface magnetic field 
of a model that can match the current observed values 
at the "real" age of 550 kyr. The model has an initial 
crustal toroidal field that reaches a maximum value of 
Stor-max(^ = 0) = 7 X 10^** G (approximately half of the 
magnetic energy is stored in the toroidal component), 
which has now decayed to Btor-max ^ 1-3 x 10^^ G. 
We have also studied the expected outb urst rate of this 
source , foll owing the same proced ure as in lPerna fc Pons! 
(|20T1 and lPons fc Perna I pOll . We found that, at the 
present stage its outburst rate is very low (w 10"'^ yr~^), 
because the magnetic field has been strongly dissipated. 

7.2. The spectral evolution during the outburst 

The spectral evolution during the outburst de- 
cay in Swift J1822. 3-1606 bears resemblance to 
that observed in other magne tars in outbu rst , no- 
tably XTEJ1810-19 7 (Halpc rn fc Gotthelf I [20051: 
Bernardini et~aI1 [20lTh . S GR0501-H45 16 (IRea et"al] 
4552 



2009D, 0X01711647-4552 (lAlbano et all 120101) and 
the "low field" SOR0418-H5729 (lEsposito et all [20T0I: 
IRea et all 120101 : iTurolla et all |201lD. The decrease in 
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Fig. 8.— Outburst model from I Pons k. Real poia l superimposed to the 1-10 keV flux decay of Swift J1822. 3-1606 (see text for details). 



flux appears, in fact, to be associated to a progressive 
spectral softening. Despite present data do not allow 
for an unambiguous spectral characterization over the 
entire outburst, evidence for a slow spectral softening 
is present in both the BB+BB and BB+PL models. In 
this respect we note that data are not consistent with 
a BB+PL fit in which the PL index is frozen to the 
same value in all observations (fitting the values of V in 
Table [3] and Figure [3] with a constant function gives a 
reduced x^ > 10). 

A BB+PL spectrum is observed at soft X-ray ener- 
gies in most magnetar sources, and is interpreted in 
terms of resonant cyclotron up-scattering of thermal sur- 
face photons by ma gnetospheric electrons in a twisted 
magnetosphere (e.g. [Thompson et al.ll2002l : iNobili et al.l 
|200 ^. In this framework, the evolution of Swift J1822.3- 
1606 is compatible with seed photons originating in a 
relatively small surface region which is heated by the 
(magnetic) event which gave rise to the outburst, mag- 
netic energy release deep in the crust (as in Lvubarskv 
I2002D and/or ohmic dissi pation of magnetospheric cur- 
rents (jBeloborodovl 120091 ). The heated region shrinks 
and cools progressively during the period covered by our 
observations (the equivalent BB radius decreased from 
~ 5 km to '^ 1 km; in the following we always assume a 
5 kpc distance) as residual heat is radiated away and the 
non-thermal component shows a progressive softening as 
the magnetosphere untwists. 

On the other hand, the spectral evolution of the source 
can be also accommodated in the framework of a BB+BB 
spectral decomposition. In this model, the thermal emis- 
sion is usually associated with two regions of different 
temperature and size which were heated during the out- 
burst. It is well possible that a single heated region is 
actually produced, but with a meridional temperature 
gradient, which can be schematized as e.g. a hotter 
cap surrounded by a warm ring, similarly t o the case 
of XTEJ1810-197 (jPerna fc Gotthelf I [20081 ) . The ab- 
sence of a non-thermal tail is not in contrast with the 
twisted magnetosphere model if the twist is small and/or 
it af fected only a limited bundle of closed field lines (see 
e.g. lEsposito et al.|[20To() , especially if the surface field 



is low, as in the present source. 

The archival ROSAT observation show that, in quies- 
cence, the source has a blackbody spectrum with kT ~ 
0.2 keV and R ~ 5 km. Although the radius is somehow 
small, it is not unreasonable to associate the ROSAT BB 
to thermal emission from the entire star surface, given 
the large errors and the uncertain distance determina- 
tion. 

If the outburst produced a heated region, which for 
concreteness we take to be a two-temperature cap, dur- 
ing the decay we witnessed a gradual shrinking of the 
hotter region (from ~ 5 km to < 1 km). The warm 
ring also shrunk and cooled down slowly during the first 
300 days after the outburst. Given the very slow spec- 
tral evolution of this component, we could obtain a good 
spectral modeling by fixing its temperature and radius 
to be constant during the first 100 days of the outburst, 
and again (at a different value) in the last 200-300 days 
(see Figure[3|) . This should be most probably interpreted 
as a gradual cooling which could not be followed in detail 
by the current observations, rather than a temperature 
jump.x 

7.3. The outburst decay and timescales 

The aggressive monitoring campaign we present here 
allowed us to study in detail the flux decay of 
Swift J1822. 3-1606, and give an estimate of its typical 
timescale. Fitting the flux evolution in the first 225 days 
after the onset of the bursting activity, we found that 
an exponential function or a power law alone cannot fit 
the data properly, since at later time (50-80 days) the 
decay slope starts to change. We found an acceptable fit 
with an analytical function of the form Flux(i) = A'l + 
ATze-^*/^) + Kst (x2=4.7/37 dof); the best values of the 
parameters are Ki = (1.76 + 0.03) x 10~^^ergs~^cm~^, 
K2 = (22.0 + 0.3) X 10-"ergs-icm-2, r = 14.6 + 0.3d, 
and K3 = (-5.2 + 0.2) x lO-^Vgs-^cm-^d-^ The out- 
burst decays of other magnetars are usually fitted by two 
components: an initial exponential or power-law compo- 
nent accounting for the very fast decrease in the first 
days or so (successfully observed only in very few cases) , 
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Fig. 9. — Perio d-period derivative diagram for all known isolated pulsars. Black dots are radio pulsars (from the ATNF Catalog; 
IManchester et an r2005\ while red symbols are all known magnetars. Asterics denote PSR J1846— 0258 and PSR 1622— 49 50, and empty 
stars are Swift J1822. 3-1606 and SGR 0418+5729. Empty grey circles are the X-ray Dim Isolated Neutron Stars (XDINS: rrurollall2D09 ;). 
The dashed line represents the value of the critical electron magnetic field. 

which is probably too last (1-10 hours) to be observable 
in most ol the cases. On the other hand, most ol the light 
curve turns out to be insensitive to the inner limit (high 
density) ol the injection region. Only the outburst tail 
(at > 50 days) is affected by this parameter, but this el- 
lect is hard to be distinguished from similar effects from 
other microphysical inputs (e.g. varying the impurity 
content ol the crust). Finally, variations ol the angu- 
lar size can be partially compensated by changes in the 
normalization lactor which at present is undetermined 
(unknown distance). This changes the volume implied 
and therefore the estimate ol the total energy injected. 
Thus we need to wait lor the lull return to quiescence, 
and combine our study with the complete analysis ol the 
pulse profile and outburst spectrum, before we can place 
better constraints on the affected volume and energetics. 



follo wed by a mu ch flat ter power-law (see IWoods et al.l 
12004 [Israel et al.ll2007t lEsposito et al.1 120081) . However. 

we note that the source has not reached the quiescent 
level yet; hence the modeling ol the outburst, and rela- 
tive timescale, might change slightly when adding lurther 
observations until the complete quiescent level is reached. 

We have also compared the observed outburst decay 
with the more physical theoretical model presented in 
Pons & Rea (2012). We have performed numerical sim- 
ulations with a 2D code designed to model the magneto- 
thermal evolution ol neutron stars. The pre-outburst pa- 
rameters are flxed by fitting the timing properties to the 
secular thermal evolution presented in section t)7.1l We 
assume that Swift J1822. 3-1606 is presently in an evo- 
lutionary state corresponding to that ol the model pre- 
sented in Figure[7]at an age ol 550 kyr. We then model 
the outburst as the sudden release ol energy in the crust, 
which is the progressively radiated away. We have run 
several ol such models varying the total injected energy 
(between lO"*" — 10^"* erg), as well as the affected volume, 
which are the two relevant parameters affecting the out- 
burst decay (coupled with the initial conditions which 
were explored in ^7.ip . The depth at which the energy is 
injected and the injection rate bear less influence on the 
late-time outburst evolution (Pons & Rea 2012). 

In Figure[5] we show our best representative model 
that reproduce the observed properties ol the decay ol 
Swift J1822. 3-1606 outburst. This model corresponds 
to an injection ol 4 x 10^^ erg cm~'^ in the outer crust, 
in the narrow layer with density between 6 x 10^ and 
6 X 10^" g cm""^, and in an angular region ol 35 degrees 
(0.6 rad) around the pole. The total injected energy was 
then 1.3 X 10'*^ erg . 

However, we must note that this solution is not unique 
and the parameter space is degenerate. Equally accept- 
able solutions can be found varying the injection energy 
in the range 1— 20x 10^^ erg cm~^ and adjusting the other 
parameters. The outer limit (low density) olthe injection 
region affects the timescale ol the rise ol the light curve. 



7.4. Radio and optical constraints 

A recent study on the emission ol radio magnetars has 
shown that all magnetars which exhibited radio pulsed 
emission, have a ratio ol quiesce nt X-ray lumino sity to 
spin-down power Lqui/irot < 1 (iRea et aT]|2012D . This 
suggests that the radio activity ol magnetars and ol ra- 
dio pulsars might be due to the same basic physical 
mechanism, while its different observational properties 
are rather related to the different topology ol the ex- 
ternal mag i ietic f ield (e.g. a dipole and a twisted field; 
lThompsonl(|2008l ). 

In the case ol Swift J1822. 3-1606, inferring the quies- 
cent (bolometric) and spin-down luminosities from our 
ROSAT data and our timing results (see §?? and ij3.2p . 

we derive Lqui/irot — 4x10^^ ergs"^/1.7x 10^° ergs~^~ 
235 . This value is in line with the^ource not showing any 
radio emission (see iRea et al.ii201 2 lor lurther details). 

Concerning the optical and infrared observations, the 
bright optical fluxes ol the sources S1-S3, much brighter 
than that ol any other SGR in outburst lor a compa- 
rable distance and interstellar extinction, as well as the 
lack ol relative flux variability, suggest that objects SI- 
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S3 are most likely unrelated to Swift J1822.3-1606. The 
hydrogen column density derived from the X-ray spec- 
tral fits (A/ji = 0.5 X 10^^ cm~^) corresponds to an inter- 
stellar extinction E(B — V ) ^ 0. 89 according to the rela- 
tion of lPredehl fc Schmi tt (1995) . Using the wavelength- 
dependent extinction coefficients of Fitz patrick (1999), 
this implies an absorption of Az ~ 1.33 and Ak ^ 0.32 
in the z and K band, respectively. Our OSIRIS upper 
limit would then correspond to an extinction-corrected 
spectral flux of < 16.3/iJy in the z band, or to an in- 
tegrated flux < 1.36 X 10~"'^'*ergcm~^ s~^. For an as- 
sumed Swift J1822. 3-1606 distance of 5kpc, this implies 
an optical luminosity L^ < 4 x 10'^"'^ erg s^"'^ during the 
outburst phase. Only very few magnetars are detected 
in the optical band, the only ones being 4U 0142-1-614, 
IE 1048.1-5937, and SG R 0501-^4516, a l l detected in 
the I band (see, e.g. iMignanil 120111 : iDhillon et"all 
[Mil). The optical flux of, e.g. IE 1048.1- 5937 during 
its recent fiaring phase (jWang et al.l 120081) corresponds 
to an /-band luminosity Lj ~ (4 ± 2) x lO^^ergs"^ 
for a n assumed distance of 3 ± 1 kpc (jGaensler et al.l 
|2005| ). Barring the difference in comparing fiuxes ob- 
tained through two slightly different filters, our optical 
luminosity upper limit is about an order of magnitude 
above the IE 1048.1—5937 luminosity. Similarly, the flux 
upper limit derived from the UKIDSS data implies for 
Swift J1822.3-1606 an IR luminosity Lr < 1.1 x lO^^ 
ergs^^ during the quiescent phase. This upper limit is 
about an order of magnitude above the computed lumi- 
nosities of the magnetars' IR counterpart in quiescence 
and, therefore, it is not very constraining. 

8. CONCLUSIONS 

We have reported on the outburst evolution of the 
new magnetar Swift J1822. 3-1606, which, despite its rel- 
atively low magnetic field {B — 2.7 x 10^^ G), is in line 
with the outbursts observed for other magnetars with 
higher dipolar magnetic fields (similar energetics, flux 
evolution and spectral softening during the decay). Fur- 
thermore, we showed that the non detection in the radio 
band is in line with its hig h X-ray conversi on efficiency 
(Lqui/irot — 235; see also iRea et al.l (|2012| ) for further 
details) . 

We studied the secular thermal evolution of 
Swift J1822. 3-1606 on the basis of the actual value 
of its period, period derivative and quiescent luminosity, 
and found that the current properties of the source can 



be reproduced if it has now an age of ^ 550 kyr, and 
it was born with a toroidal crustal field of 7 x 10"'^'* 
G, which has by now decayed by less than an order of 
magnitude. 

The position of Swift J1822.3-1606 in the P~P diagram 
(see Figure|9]) is clos e to that of the "low" field magne- 
tar SGR 0418-1-5729 (|Rea et al.ll2010t ). Although the fact 
that both have a sub-critical dipole field is not relevant 
per se, and the dipolar field in Swift J1822. 3-1606 is at 
least four times higher than SGR 04184-5729, it is worth 
to stress that the discovery of a second magnetar-like 
source with a magnetic field in the radio-pulsar range 
strengthens the idea that magnetar-like behavior may 
be much more widespread than what believed in the 
past, and that it is related to the intensity and topology 
of the internal and surface toroidal cor nponents, rather 
than only to the su r face dipolar field ( Rea et aJ] 120101 : 
IPerna fc PoiisllmllTilrolla et al.llMlL 

Monitoring the source until its complete return to qui- 
escence will be crucial to disentangle: 1) its complete 
spectral evolution during the outburst decay, 2) the pos- 
sible presence of a second derivative of the rotational 
period, possibly due to the source timing noise, 3) re- 
fine further the modeling of the outburst and the surface 
region affected by this eruptive event. 
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